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Abstract 

Maxwell's equations for the time-dependent electromagnetic field 
in a homogeneous chiral medium are reduced to a single quaternionic 
equation. Its fundamental solution satisfying the causality principle is 
obtained which allows us to solve the time-dependent chiral Maxwell 
system with sources. 

1 Introduction 

We consider Maxwell's equations for the time- dependent electromagnetic 
field in a homogeneous chiral medium and show their equivalence to a sin- 
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gle quaternionic equation. This result generalizes the well known (see [T3] . 
0; 501) quaternionic reformulation of the Maxwell equations for non-chiral 
media. Nevertheless the new quaternionic differential operator is essentially 
different from the quaternionic operator corresponding to the non-chiral case. 
We obtain a fundamental solution of the new operator in explicit form sat- 
isfying the causality principle. Its convolution with a quaternionic function 
representing sources of the electromagnetic field gives us a solution of the 
inhomogeneous Maxwell system in a whole space. 

2 Maxwell's equations for chiral media 

Consider time-dependent Maxwell's equations 

rot 'E{t, x) = -Si's it, x), (1) 

iotli{t,x) = dt'D{t,x) + y{t,x), (2) 

div 'E{t, x) = div H{t, x) = (3) 

with the Drude-Born-Fedorov constitutive relations corresponding to the chi- 
ral media 0, [H], US 

'B{t,x) = fi(H{t,x) + 13 lot H{t,x)), (4) 

'D{t,x) = e(E{t,x) + PTot^{t,x)), (5) 

where f3 is the chirality measure of the medium. f3, e, are real scalars 
assumed to be constants. Note that the charge density p and the current 
density j are related by the continuity equation dtp + div j = 0. 

Incorporating the constitutive relations (jl)), (0) into the system (P)-© we 
arrive at the main object of our study, the time-dependent Maxwell system 
for a homogeneous chiral medium 

rot H{t,x) = e{dt'E{t,x) + pdtTot^{t,x)) + y{t,x), (6) 
rot E (t, x) = -p{dtHit, x) + (3dt rot ^(t, x)), (7) 
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divE(t,x) = div H{t,x) = 0. (8) 

Application of rot to (jHI) and ((Zj) allows us to separate the equations for 
E and H and to obtain in this way the wave equations for a chiral medium 

(9) 



rot rot H + efid^H + 2(3eixdf rot H + P'^efidf rot rot H = rot j . (10) 

It should be noted that when /3 = 0, and (fTUjl reduce to the wave 
equations for non-chiral media but in general to the difference of the usual 
non-chiral wave equations their chiral generalizations represent equations of 
fourth order. 



3 Some notations from quaternionic analysis 

We will consider biquaternion- valued functions defined in some domain Q C 
M'^. On the set of continuously differentiable such functions the well known 
Moisil-Teodoresco operator is defined by the expression D = + + 
(see, e.g., 0), where ik, k = 1,2,3 are basic quaternionic imaginary 
units. Denote = D + a, where a G C and Ima > 0. The fundamental 
solution for this operator is known j9| (see also pil] ): 



}Ca{x) = — grad 0a(x) + aQa{x) = (a H 5- — zct-j— j-)0Q(x), (11) 

\x\ \x\ 

where i is the usual complex imaginary unit commuting with ik, x = Yl'k=i -^^ik 
and 0a (x) = — ^'^1^1 . Note that /Cq, fulfills the following radiation condition 
at infinity uniformly in all directions 

XX 1 
(1 + 1 — i) ' ^a{x) = o{—) , when |x| — > 00 (12) 

\x\ \x\ 

which is in agreement with the Silver-Miiller radiation conditions [Hj. 
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4 Field equations in quaternionic form 

In this section we rewrite the field equations from Section |21 in quaternionic 
form. 

Let us introduce the following quaternionic operator 

M = py/eJldtD + y/ej2dt - iD (13) 
and consider the purely vectorial biquaternionic function 

V{t,x) = E{t,x) -i^H{t,x). (14) 

Proposition 1 The quaternionic equation 

MV{t, x) = -^f 7(t, x) - /3^d,pit, x) + '-P^ (15) 

is equivalent to the Maxwell system (0)-(0), the vectors E and H are solu- 
tions of (0)-(0) if and only if the purely vectorial biquaternionic function V 
defined by [T^ is a solution of / f73|) . 

Proof. The scalar and the vector parts of (jlSp have the form 
-(3^dtdiv^+J^diYH+i{dw^+PtidtdiYH) = -pJ^dtp+j, (16) 



/3v^5trot^+^/£//at^-Y'^rot^-z(rot^+/?/iatrot^+//<9t^) = -y ~7 

(17) 

The real part of ()17p coincides with © and the imaginary part coincides 
with (|7j). Applying divergence to the equation (fTTj) and using the continuity 
equation gives us 

dt div H = and dt div E = -dtp. 

Taking into account these two equalities we obtain from ()l(j|l that the vectors 
E and H satisfy equations (jHI). ■ 
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It should be noted that for (3 = from (|13|) we obtain the operator which 
was studied in [7j with the aid of the factorization of the wave operator for 
non-chiral media 

sjjd'f — Ax = {y/ej2dt + iD){y/eJldt — iD). 

In the case under consideration we obtain a similar result. Let us denote by 
M* the complex conjugate operator of M: 

M* = py/eJIdtD + y/ejldt + iD. 

For simplicity we consider now a sourceless situation. In this case the equa- 
tions (jni) and (fTUj) are homogeneous and can be represented as follows 

MM*U{t,x) = 0, 

where U stands for E or for H . 



5 Fundamental solution of the operator M 

We will construct a fundamental solution of the operator M using the re- 
sults of the previous section and well known facts from quaternionic analysis. 
Consider the equation 

{(3^/ejldtD + y/ejldt — iD)f{t, x) = S{t, x). 

Applying the Fourier transform JF with respect to the time-variable t we 
obtain 

{Py/ejliujD + y/ejliu! — iD)F{uj, x) = 6{x), 

where F{u,x) = J-'{f(t,x)} = fit,x)e^^'^'^dt. The last equation can be 
rewritten as follows 

{D + a){(3y/ejluj — l)iF{uj, x) = 5{x), 

where a = j^=^zi - The fundamental solution of is given by (fTT|) . so we 
have ^ ^ 

{(iJeJiu — l)iF(u}, x) = (a -\ ^ — ia— — )9q(x), 

\x\ \x\ 

from where 
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F{lo,x) 



IWEflUJ 



{(3^u - 1) 



IX \ IX 

+ 

\x\ 



\x\ Py/eJIuj - 1 



i\x\ 



An \x\ 



We write it in a more convenient form 

1 



F{uj,x) 



{u — ay 



jA(x) 



I \ ic(x) 

B(x) E(x)e—. 



bj — a 



where a = c (x) = E (x) 



i\x\ 

4:ir\x\ ' 



1 / ix\ X 



In order to obtain the fundamental solution f{t, x) we should apply the in- 
verse Fourier transform to F{uj,x). Among different regularizations of the 
resulting integral we should choose the one leading to a fundamental solution 
satisfying the causality principle, that is vanishing for t < 0. Such an election 
is done by introducing of a small parameter y > in the following way 



f{t,x) = \imj^-^ {F{z,x)} 



(18) 



where z = u — iy. This regularization is in agreement with the condition 
Ima > 0. We have 



:f-'{f{z,x)} = ^ r f-^ 



rA (x) + — - — B (x) ] E (x) e"-"^ e'^^'du 



where ay = a + iy. Expression (fTUI) includes two integrals of the form 



(19) 



27r (cu - ay)' 



:du, k = 1,2 



where c = c(x) . We have 



2n 



E 

j=0 



Uc 



J roo iujt 



3- J-oo {oJ — ay) 



j+k 



(20) 
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Denote 

h,j{t) — / - — 
J-oo [uj — ayY 

For k = 1 and j = we obtain (see, e.g., ^ Sect. 8.7]) 

hflit) = 27riH{t)e'"'y 

where H is the Heaviside function. For all other cases, that is for k = 1 
and j = 1, oo and for k = 2 and j = 0, oo we have that j + k > 2 and the 
integrand in ()20j) has a pole at the point ay of order j + k. Using a result 
from the residue theory |3J Sect. 4.3] we obtain 

Ik j{t) = 2711 Res for t > and j + k>2. 

Consider 

Res -r = - i- hm , ^ e--*= V ^ , „ for t > 

(a;-aj,y+'' (j + A; - 1)! ..-a, ^cj^+fc-i + A; - 1)! 

and j + k>2. 

For t < we have that Ik,j{t) is equal to the sum of residues with re- 
spect to singularities in the lower half-plane y < which is zero because the 
integrand is analytic there. Thus we obtain 



4,(t)=2vrz//(t)^A_L__e-.* 
Substitution of this result into ()20|) gives us 



j=o •' ■' j=0 -J ^-^ ' > 

Now using the series representations of the Bessel functions Jq and J\ (see 
e.g. jUI Chapter 5]) we obtain 

h = iH{t)e'''y'Jo (2^/a^ and h = -H{t)^e'''^'Ji (2^0^ . 

Substituting these expressions in p9|) and then in (fTHj) we arrive at the 
following expression for /: 
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Finally we rewrite the obtained fundamental solution of the operator M in 
explicit form: 

f{t,x) = Hit)^ (jCMJo 

Let us notice that / fulfills the causality principle requirement which 
guarantees that its convolution with the function from the right-hand side of 
(fT^ gives us the unique physically meaningful solution of the inhomogeneous 
Maxwell system (0)-® in a whole space. 
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